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Abstract 


‘A  geometric  programming  model  is  proposed  to  determine  the 
optimal  design  of  the  CPU  and  its  matching  storage  hierarchy.  The 
objective  function  is  the  maximization  of  system  reliability  subject 
to  performance  and  budgetary  limitations.  Examples  illustrating  the 
use  of  the  model  are  presented. a 
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1.  Introduction 


An  important  area  of  concern  during  the  design  of  a  computer 
system  is  the  design  of  its  storage  subsystem.  Several  optimization 
models  for  the  design  of  storage  hierarchies  are  available  [3,4,8,12]. 
These  models  typically  optimize  system  performance,  measured  by 
throughput  or  average  hierarchy  access  time,  subject  to  budgetary  lim¬ 
itations.  In  the  previous  efforts,  no  consideration  is  given  to  the 
reliability  issues.  In  systems  designed  for  avionics,  space,  and  cer¬ 
tain  military  applications,  reliability  is  of  utmost  concern.  In  this 
paper,  we  develop  an  optimization  model  for  linear  storage  hierarchies 
with  the  objective  of  maximizing  system  reliability  subject  to  a  cost 
and  a  performance  constraint.  The  capacities  of  various  memory  levels 
are  the  decision  variables.  An  extended  model  also  includes  the  CPU 
speed  as  a  decision  variable. 

Chow  [3,4],  considered  the  design  of  a  linear  storage  hierar¬ 
chy  with  the  objective  of  minimizing  the  average  access  time  subject 
to  a  cost  constraint.  The  present  paper  is  an  extension  of  Chow's 
model  to  include  reliability  considerations.  Other  related  work  is  by 
MacDonald  and  Sigworth  [8] ,  and  Trivedi  and  Sigmon  [12] . 
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2 .  Model  Development 

The  storage  hierarchy  consists  of  n  levels,  •••» 

as  shown  in  Figure  1.  The  hierarchy  management  follows  the  staging 
rules  proposed  by  Gecsei  [6] .  The  execution  is  out  of  the  fastest 
level  M,  which  may  be  thought  of  ns  a  cache.  Information  transfers 
are  between  adjacent  levels  only.  By  convention,  the  higher  the  level 

(smaller  the  index)  in  the  hierarchy,  the  smaller  its  capacity.  Rules 

of  operation  of  such  a  hierarchy  are  described  in  16,10). 

The  capacity  of  memory  level  M^,  will  be  denoted  by  (bits), 
and  the  block  size  by  s^  (bits) .  Since  the  entire  address  space  of 
the  program  must  be  contained  in  level  n,  we  assume  that  the  capacity 
of  that  level,  K^,  is  a  fixed  input  parameter  to  the  model.  The 
decision  variables  are  the  capacities  K^,  K2t  •••»  *^n-l* 

The  behavior  of  the  program  can  be  characterized  by  a 
reference  string.  We  assume,  however,  that  the  program  behavior  and 
the  effect  of  storage  management  strategy  are  compacted  into  a  single 
function  called  the  success  function  which  gives  the  probability  that 
a  storage  reference  from  the  CPU  is  found  in  a  given  level.  The  suc¬ 
cess  function  for  level  i  depends  on  the  capacity  of  level  i  (K^^)  / 

the  block  size  of  level  i  and  all  of  the  lower  levels 

(Si,  •••»  Sj^)  ,  and  the  block  replacement  algorithm.  If  the 

hierarchy  management  rules  satisfy  certain  reasonable  properties  [6], 
it  can  be  shown  that  the  success  function,  H^,  depends  only  upon  the 
capacity  and  the  block  size  s^.  That  is, 

Hj^  =  Hj^(Kj^;s^)  1  =  1,2,. ..,n  . 
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The  semicolon  is  used  above  to  emphasize  the  fact  that  is  a 
decision  variable  while  s^^  is  a  fixed  parameter.  The  miss  ratio  func¬ 
tion  is  defined  to  be 


F.  =  1  -  H. 

1  1 

Let  t^  be  the  average  time  to  transfer  a  block  of  size 
from  the  i^*^  level  to  level  i-1.  The  access  times  t^,  t^^  are 

considered  fixed  parameters  in  our  model.  Due  to  the  linear  organiza¬ 
tion,  the  time  to  fetch  a  block  from  level  i  and  percolate  it  (or 
appropriate  sub-blocks)  up  to  level  1,  denoted  by  T^,  is  given  by 


The  average  access  time,  T,  of  the  entire  hierarchy  is  the 
weighted  average  of  the  access  time  to  various  levels,  that  is, 
n 

T  =  Z  h.T. 
i=l  ^  ^ 

The  weight  h^  is  the  probabil^ity  of  a  hit  to  level  i  and  misses  to  all 
the  previous  levels  j  <  i.  Thus 


h.  =  H.  -  H.  .  =  F.  ,  -  F. 
1  1  1-1  1-1  1 


1—1,2, . . .n 


and 


Therefore,  the  expression  for  T  can 


=  "i  +  . 

1=2 

n 

=  ’'l  ^  .^.’^i^i-l^^i-l'^i-l’ 
1  =  2 


by  convention, 
be  simplified  to 
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The  cost  of  the  memory  level  is  modeled  by  a  posynomial 
function  of  its  capacity,  that  is 


DEVCOST.  =  5-  c.  -K  . 

1  j  ''ID  1 


j  2.  f  , 2, . . .  ,n 


The  total  svstem  cost  is  then  given  by  Z  DEVCOST.. 

i=l  ^ 

The  failure  rate,  of  memory  level  i  is  a  posynomial  func¬ 

tion  of  its  capacity,  that  is 

B.  . 

^i  ?  '^ij^i  ^  '  ^ij  ^  ®  '  i=l»2,...n  . 

Note  that  the  military  standard  M1L-217B  [9]  suggests  that  is  a 

linear  function  of  capacity  K^.  A  failure  in  any  level  is  assumed  to 
cause  the  entire  hierarchy  to  fail.  Thus,  from  the  reliability  point 
of  view,  it  is  a  series  system.  If  we  further  assume  that  the  failure 
rate  of  each  level  is  independent  of  its  age  (or  equivalently,  the 
lifetimes  are  exponentially  distributed),  the  hierarchy  failure  rate, 
X,  is  also  time-independent,  and  is  given  by  [1], 
n 

X  =  S  X.  . 

i=l  ^ 

For  such  a  system,  minimizing  the  failure  rate  is  equivalent  to  maxim¬ 
izing  the  reliability  [1]. 

We  are  now  ready  to  define  the  optimization  problem  for  the 
design  of  the  storage  hierarchy.  Assume  that  the  hierarchy  access 
time  is  not  to  exceed  T«  and  the  system  cost  is  constrained  by  BUDGET. 


min  X  =  z  2  d.^K.- 
i=l  j  ^ 


such  that 
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n  . 

Z  Z  c  -  -K.  -*  <  BUDGET 
i=l  j  ^ 


K^  >  0  ,  i=l,2, . . . ,n-l 


(lb) 

(Ic) 

(Id) 


If  we  assume  that  the  miss  ratio  function,  F^,  is  a  posynomial 

function  of  capacity  K.  (that  is, 

■^ii  " 

=  Z  ■’  ,  '  1=1, 2, .  .  .  ,n-l)  ,  then  the  design  problem 

j  J  3 

(1)  above  is  a  standard  geometric  programming  problem.  From  the 
theory  of  geometric  programming  [13] ,  we  get  the  following  result: 


Theorem  _1: 

Any  relative  minimum  of  the  design  problem  (1)  is  also  its  global 
minimum. 


The  implication  of  this  result  is  that  any  standard  technique 
for  nonlinear  optimization  will  locate  the  global  minimum  of  the 
design  problem.  Alternatively,  geometric  programming  techniques  can 
be  used  to  solve  the  problem.  In  the  next  section,  we  consider  a  sim¬ 
plification  of  this  problem  so  that  a  closed  form  solution  can  be 
obtained . 


3.  Simplified  Design  Problem 

In  the  previous  section,  the  miss  ratio  function,  the  device 
cost  function,  and  the  failure  rate  function  were  assumed  to  be  quite 
general  posynomial  functions.  In  this  section  we  make  the  following 
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simplifying  assumptions  (that  is,  restrict  the  number  of  terms  in  the 
posynomials) : 


f  i  — 1/2/.**/!!  1 


^i  ■  '^Oi  ^li^i  '  i=l»2,...,n 


/i 


DEVCOST.  =  c„.  +  C,.K.  ,  i=l,2, 

1  Oi  li  1 


/n 


(one  term  only) 
(two  terms  only) 
(two  terms  only) 


Note  that,  since  K  is  fixed,  X  and  DEVCOST  are  fixed.  Define 

n  n  n 

n  B 

d.  =  S  d-.  +  d,  k'  "  , 

0  0i  In  n 


Now  the  design  problem  (1)  can  be  rewritten  as 

n-1  p. 

min.  d-  +  Z  d,.K^ 

0  i  =  l  ^ 


(2a) 


n 


-^i-1 


1  .^2  ^  ~  ° 


n-1 


<  BUDGET  -  Cq 

K.  >  0  , 

•  •  •  f  rj”  ^ 

(2b) 

(2c) 

(2d) 


With  these  assumptions,  it  can  be  shown  that  [11]  the  solution 
to  design  problem  (1)  can  be  decomposed  into  two  steps.  The  first 
step  is  to  obtain  the  solution  to  the  following  subproblem; 


n-1 


min.  dp  +  £ 

i=l 


Pi 


(3a) 


ti  +  .J^^iSi-iKi-l  <  I’d 


(3b) 
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>  0  ,  i-1 , 2 , . . . ,n-l  (3c) 

In  particular,  we  have  remover!  the  cost  costraint.  Assume  that  a 
solution  to  problem  (3)  is  given  by  K^,  K2t 
n-1  /. 

i=l  ^  -  0 

then  this  is  also  the  solution  to  the  overall  design  problem  (2).  On 
the  other  hand,  if 
n-1 

Z  C,  .K .  >  BUDGET  -  c^, 

i=l  ^  ^ 

then  the  original  design  problem  (2)  possesses  no  feasible  solution. 
In  other  words,  the  requirements  of  the  design  problem  are  incompati¬ 
ble.  Due  to  this  result,  we  will,  henceforth,  omit  the  cost  con¬ 
straint  from  the  problem  specification.  We  will  assume  that  the 
requirements  are  compatible. 


With  the  further  restriction  that  c(^  =  c(  and  =  p  for  all  i, 
the  following  closed  form  solution  to  design  problem  (3)  can  be 
obtained  [ 11] . 


As  an  example,  we  take  the  input  parameters  shown  in  Table  1 
for  the  design  of  a  3-level  storage  hierarchy.  The  data  on  the  miss 
ratio  function  is  derived  from  [8].  All  other  input  parameters  are 
based  on  the  data  reported  in  [2,5].  The  resulting  optimal  capacities 
are  also  shown  in  Table  1. 
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Next  we  consider  the  problem  with  the  same  parameters,  except 
that  the  maximum  allowable  hierarchy  access  time,  is  to  be  varied. 
The  result  of  such  a  sensitivity  analysis  is  shown  in  Figure  2.  We 
have  plotted  the  optimal  hierarchy  MTTF  (equal  to  the  reciprocal  of 
the  system  failure  rate)  as  a  function  of  the  maximum  allowable 
hierarchy  access  time  T^.  In  Figure  3,  we  have  plotted  the  cost  of 
the  optimally  designed  hierarchy  as  a  function  of  T^. 

4 .  Including  the  Selection  of  CPU  Speed 

The  model  of  linear  storage  hierarchy  discussed  up  to  this 
point  does  not  address  the  selection  of  CPU  speed.  We  have  also 
ignored  the  effect  of  CPU  instruction  execution  delay  from  performance 
considerations.  Let  tp  be  the  average  instruction  interpretation  and 
execution  time  (i.e.,  1/tQ  is  the  instruction  execution  rate  assuming 
an  infinitely  fast  memory  hierarchy) .  Let  A  be  the  average  number  of 
memory  references  per  instruction.  Gecsei  and  Lukes  [7]  estimate  that 
A  =  2  while  Snow  and  Siewiorek  [2]  estimate  that  A  =  1.162. 


The  average  execution  time  per  instruction  is  now  given  by 
to  t  ''Itl  +  j2tiF(K._j)) 


Apart  from  the  memory  capacities  K^,K2f • . . we  consider 
tp  also  to  be  chosen  by  optimization.  We  assume  that  tp  is  a  power 
function  of  the  complexity,  Gq,  of  the  CPU: 
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Gq  may  be  equated  to  the  gate  count,  or  any  related  measure  of  com¬ 
plexity.  The  CPU  failure  rate,  Xq»  is  also  a  power  function  of  the 
complexity: 


^0  *^00  '^lO^O 


Po 


Finally,  the  CPU  cost  is  modeled  by 


‘^OO  ^10^0 


The  design  problem  (1)  can  be  redefined  as  follov;s: 


Pi  Po 

min.  d  +  d^.K.  +  d^gG 

1  =  1 


(5a) 


®10=0  ^  i  TES  -  «00 


y  n  y 

=00  *  =10°0°  *  .!/=0i  * 


(5c) 


Gq  >  0 

K .  >  0 
1 


(5d) 

(5e) 


As  in  section  3,  we  use  the  miss  ratio  function 
"'^i 

=  a^K^  ,  i=l,2, . . .  ,n-l 

lER  is  the  minimum  required  instruction  execution  rate  for  the  system 


consisting  of  the  CPU  and  the  storage  hierarchy.  Also, 


Pn 


d_  =  Z  d_.  +  d,  K 
0  Oi  In  n 


As  before,  we  can  solve  the  subproblem  ignoring 


the  cost  constraint.  Assuming  that  c(^  =  c(,  =  p  for  i-1, 2, . .  .  ,  n-1 , 

it  can  be  shown  that  the  CPU  delay  tp  must  satisfy  the  equation  [11]: 


where 


p  (c(+p)C(Q 


This  is  a  nonlinear  equation  with  one  unknown  and  can  be  solved  using 
the  standard  iterative  techniques  to  obtain  the  optimal  value  of  t^. 
Now  the  optimal  values  of  the  memory  capacity  ( i  =  l , 2 , . . . , n-1)  can 
be  obtained  from  equation  (4)  after  substituting 


A  4eR 


- 


for  the  term  Tq  -  tj^. 


As  an  example,  assume  we  want  to  design  a  system  which  is 
capable  of  executing  at  0.2  MIPS.  Other  input  parameters  are  shown  in 
Table  2.  The  parameter  values  are  based  on  data  reported  in  [2,5,8]. 
The  optimal  solution  is  also  shown  in  Table  2.  The  design  roughly 
corresponds  to  a  typical  PDP-11/60  configuration  [2] . 


Next,  we  perform  sensitivity  analysis  varying  the  specified 
value  of  the  lER.  The  system  MTTF  and  the  system  cost  are  plotted  as 
functions  of  lER  in  Figures  4  and  5,  respectively. 


5 .  Concluding  Remarks 

We  have  developed  a  reliability-oriented  design  model  for 
linear  storage  hierarchies.  The  decision  variables  are  the  capacities 
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of  each  memory  level  and  the  speed  of  the  CPU.  The  design  problem  is 
set  up  as  a  geometric  programming  problem  with  the  objective  of  maxim¬ 
izing  system  reliability  subject  to  a  cost  and  a  performance  con¬ 
straint.  Any  relative  optimum  is  a  globally  optimum  solution  to  the 
design  problem.  Closed  form  expressions  for  the  device  capacaties  are 
obtained.  Examples  illustrating  the  usefulness  of  this  model  are  also 
presented . 
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Table  1  :  Example  of  a  Linear  Storage  Hierarchy  Design 

INPUT  PARAMETERS  : 

Number  of  Levels,  n  =  3 

Desired  Hierarchy  Access  Time,  T^  =  2.41  psec 

Address  Space  Size,  =  500K  words 

Miss  Ratio  Function,  F.  =  a.K.^ 

1  11 

where  ,  a^^  =  576 

,  32  =  7056 

,  c(  =  1.449 


Memory  Level 

Access  Time 
per  word 

Block  Size 
(words) 

t . 

1 

d,  . 
li 

per  lO^hr/word 

^li 

$/word 

1  (cache) 

1  0.1  ps 

1  32  ; 

1  0.1 

1  0.066 

1  1.60 

2  (main  mem) 

1  1  ps  1 

[  128  1 

I  32  1 

[  0.004625  1 

1  0.160 

3  (back  up  ) 

(  memory  ) 

10  ps 

- 

1280 

- 

0.016 

1 

dg  =  0  ,  Cq^  =  0  ,  dp^  =  0  ,  p.  =  1  ,  =  1  for  all  i  . 


OUTPUT  (OPTIMAL  DESIGN) 


Memory  Level 


Optimal 

Contribution 

Contribution 

Capacity 

to  Access 

to  Failure 

K .  (words) 

Time 

Rate  , 

X 

(ps) 

per  10°hr 

1.2K 


0.1 


79.2 


Contribution 
to  cost 

($) 

1920 


2 

1  44.35K  1 

0.641  1 

1  204.7  1 

7080 

3 

1  500K  1 

1.669  1 

1  -  1 

MCO 

Totals 


2.410 


283.9 


17,000 


Table  2  : 


Example  Design  of  a  CPU  and  its  Matching  Storage  Hierarchy 

INPUT  PARAMETERS  ; 

Number  of  levels,  n  =  3 

Desired  Instruction  Execution  Rate,  lER  =  0.2  iMIPS 

Address  Space  Size,  K  =  500  K  words 

n 

Miss  Ratio  Function,  F.  =  a.KT*^ 

1  11 

where  ,  a^^  =  576 

,  32  =  7056 

,  c(  =  1.449 

Number  of  Storage  References  per  Instruction,  A  =  1.162 
Bq,o  =  0  ,  c(q  =  1  ,  0  ~  2.298  x  10^  gate-psec 

d^Q  =  0  ,  Pq  =  1  ,  =  0.00988  failures  per  gate  per  10®  hours 

Coo  $20,000  ,  /g  =  1  ,  CjQ  =  0.61  $  per  gate 

All  other  parameters  are  equal  to  their  corresponding 
values  in  Table  1. 

OUTPUT  (OPTIMAL  DESIGN) 


Memory  Level 

Delay, t^ 

(ps) 

Complexity 

Contribution 
to  Failure 
Rate  c 

per  lO^hr 

Contribution 
to  Cose 

($) 

0 

(CPU) 

1.915 

11971 

Gate 

Count, 

^0 

118 

27,300 

1 

(cache) 

0.1 

l.llK 

72 

1780 

2 

(main  mem) 

32 

41. 4K 

Capac . 

191 

.  ^  ^  ^  ^  1 

6620 

3 

(backup 

memory) 

1280 

500K 

(words) 

8000 

_ 

Totals 

381  1 

43,700 

i/LflMBDR  xlO 


idt*  lij  iBoi 


FIGURE  3,  Hierari 
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3.21 


S.O'l 


FIGURE  5 


5.62  6.83 

lER  CMIPS)  xlO- 


System  Cost  vs.  Instruction  Execution  Rate 


9.25 


